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ABSTRACT 

For the first time, we have used RHESSI's spatial and energy resolution to 
determine the combined chromospheric and coronal density profile of the flaring 
solar atmosphere in a statistical manner, using a dataset of 838 flares observ- 
able in hard X-rays above 25 keV. Assuming the thick-target beam model, our 
"average flaring atmosphere" was found to have density scale heights of 131±16 
km at low altitudes (chromosphere, up to ps1-1.5 Mm above photosphere), and 
of 5-6 Mm at high altitudes (corona, above ^2-3 Mm). Assuming a unit step 
change in ionization level, modeling yields a height of 1.3±0.2 Mm for the tran- 
sition between fully neutral to fully-ionized atmosphere. Furthermore, centroids 
of emission above 50 keV, produced by electrons of similar or higher energies, 
are located mostly in a small region ~0.5 Mm in vertical extent, where neutral 
densities are beyond 3xl0 13 cm -3 . 

Subject headings: Sun: flares - Sun: particle emission - Sun: X-rays - Sun: Corona 
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Introduction 



The density structure of the Sun's atmosphere in the vicinity of the transition region is 



(see e.g. 



known, at least 


for f 


Vernazza et al. 


1981 



Fontenla et al. 



1993 



Gabriel 



1976 



Ewell et al 



19931 ) 



More r e cently lAschwanden et all (120021 ); iLiu et all (120061 ): Kontar et all (120081 ): 



Prato et al. 



( 120091 ). have attempted to derive the chromospheric density structure with the 



use of hard X-ray (HXR) emi ssion from f 



the thick-target beam model ( IBrown 



are footpoints at di fferent energies, assuming 



Aschwanden et al. 



(I2OO2J) 



19711 ). To fit their data 
have ass umed that th e dens ity has a power-law shape with altitude above photosphere 
whereas 



Kontar et al. 



( 120081 ). working at slightly higher energies (and hence, dee per in t he 



chromosphere), have assu med an exponential shape. Uniform target ionization ( IBrown 



1973 



Kontar et al. 



20021 ) was assumed in both case (fully-ionized for the for mer, fully 



Ewell et al. 



neutr al for the latter). The Caltech Irreference Chromospheric Model (CICM, 
19931 ) supports a two-exponential atmosph ere, with the l ower c omponent's scale height 



closely corresponding to th e one derived by 



( 11981 



thereafter VAL) and Fontenla et al. 



Kontar et al. 



( I2008h . and to the 



Vernazza et al. 



1993L thereafter FAL) models (i.e. ~130 km). 



These previous studies used single events. 



Matsushita et al. 



ril992j ) and ISatol (120061 ) 



have statistically derived the altitude differ ence between 



at different energies using Yohkoh HXT's (IKosugi et al. 



spatia lly-averaged HXR emissions 



19921 ) four channels, covering 



an energy range between 14 to 93 keV. We used the same technique of determining 
spatially-averaged centroids (center of mass) at different energies, and in a similar energy 
range, but with flares obs erved by the Ramaty High Energy Solar Spectroscopic Imager 
(RHESSI, Lin et al. 



20021 ). which has a much higher spectral resolution (~1 keV). This will 
allow us to carry the data analysis one step further, deducing densities from the regions 
where non-thermal HXR emission is observed, assuming a thick-target beam model for 
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bremsstrahlung emission. 

We will first derive densities using a simple, direct method, and then try to fit a 
double-exponential with unit step ionization change density model. 



2. Observations 

2.1. Data selection 

We have selected all flares in the RHESSI flare listQ, between the start of the RHESSI 
mission on 2002 February 5 and 2010 January 1. Accompanying the flare list are "Quicklook 
Images" in different energy bands (3-6, 6-12, 12-25, 25-50, 50-100, and 100-300 keV), 
typically accumulated over two minutes. Quicklook Images in a certain energy band are 
made only if they were deemed to "reliably image" the X-ray source, that is, that an 
X-ray source was observed to be roughly at the same location in a majority of RHESSI 
sub-collimators. We examined all events for which such images above 25 keV existed. There 
were 838 of them. 



2.2. Derivation of altitude-energy relationship 



We assume that non-thermal electrons precipitate along magnetic field lines that 
radially extend from the photosphere to the corona. They propagate from an acceleration 
region somewhere in the co rona, and lose ener gy and emit bremsstrahlung HXR as per 



the thick-target model (e.g. 



Brown et al. 



20021 ). Under these assumptions, the difference 
in source altitude, Ah i: between energy e$ and a reference energy e rP f can be used t o 



approximate the average density between these two points (see e.g. 



Brown et al. 



2002 



, and 



^■http:/ /hesperia. nasa.gov/hessidata/dbase/hessi_flareJist.txt 
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Appendix [B]): The idea is that for an accelerated electron distribution with a negative 
power-law index propagating towards region of higher densities, non-thermal emission at 

2 

energy e spatially peaks where electrons have crossed a column density N « with K a 
constant. Hence, knowing the distance s between peaks at emission e\ and e^-, one can get 



N 2 -N! 
s 



an average density between the two peaks of emission: n 
iV = Kj7 is on ly a handy approximation: as discussed in iBrown et al. 



n actuality, 



Xu et al. 



(J2002|) (see also 



20081 ) . there is a dependence on the spectral distribution of injected electrons, and 
also on the density profile of the medium through which the electron beam is propagating. 
While both dependences can easily be taken into account when using forward-fitting 
techniques on singular events, both vary from flare to flare, and we have thus decided to 

2 

use the N = 4^ relationship throughout our statistical study. 



2002|), accumulated over 



We have computed the RHESSI visibilities (iHurford et al. 
three minutes around peak HXR flux, in the following energy bands: 6-10, 10-15, 15-20, 
20-25, 25-30, 30-35, 35-40, 40-50, 50-60, 60-70, 70-80, 80-90, 90-100. We have used 
the hsi_vis_fwdfit.pro routine to find positions of centroid at different energies, using 
subcollimators (SCs) 3-9. This method basically yields the position of the centroid of the 
flux (i.e. its center of gravity). 

The software also yields error bars on centroid positions. Typically, higher energies, 
having less count statistics, yield larger error bars. We omitted cases where the software 
did not converge to a solution. 

Once the centroid positions are determined at all energy bands and for all flares, we 

1 ene r gies E j and 



Satol (J2006|) (and 



statistically determine the altitude difference, A/ij, between emission at a 
a reference energy e re f. The A/ij are derived using the same method as in 
described in more mathematical details in Appendix . Figure [T] shows plots of R vs. AR, 
with R the projected (on the usual plane of the Sun, perpendicular to the observer's line 
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of sight) distance between HXR emission and Sun center, and AR the projected distance 
between the centroid of emissions at different energies. The slope ^ is obtained by linear 
fitting with a zero intercept coordinate (i.e. going through the origin). Because high energies 
have lesser statistics, leading to poorer positional accuracy of the centroid, and because 
there is no potential contamination in the determination of non-thermal centroid position 
by thermal emission below 35 keV, (although other effects exist and are discussed in 
Section [273]) . the reference energy was chosen to be 35 keV throughout. Note that there are 
many outliers, but that their error bars are generally larger, thus have a lesser contribution 
to the determination of Finally Ahi=r s ^, with r s =solar radius (see Appendix IA~]) . 

From the Ahi values obtained in the preceding paragraph, the plot in Figure |2] (top) 
is constructed. It shows the altitude differences between emission at reference energy 
e re /=35 keV and emission at other energies (we have also looked at other choices of e re f, 
and found similar results, albeit with larger error bars when e re f was higher than 35 keV). 
Also shown are curves derived from subsets of our data: blue: using only events with 
maximum pileup error <20%, and magenta: using only events with maximum pileup error 
>20 %. Pileup issues will be discussed further on. The results presented in Figure [2] (top) 
show a systematic decrease in altitude as a function of energy, despite somewhat largish 
error bars at high energies. Notice at low energies the departure from a strictly monotonic 
relationship. We attribute this to contamination of our sample by thermal emission at low 
energies. 

We have checked the reliability of results in Figure [2] (top) by varying certain selection 
criteria and redoing the procedure. We have tried to select events with quicklook images 
above 50 keV (as opposed to 25 keV), tried other subcollimator combinations (e.g. 1-9), 
tried varying accumulation time intervals, different energy bands, and different e re f'- the 
curve essentially remains the same, though generally with larger error bars. We have also 
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Fig. 1. — Abscissa: R the projected distance between Sun center and flare emission, ordi- 
nates: AR, the projected distance between emission centroid at 35-40 keV and at the energy 
band indicated on top of each plot, with error bars. 



verified that a histogram of the normalized residuals from the ^ linear fitting are normally 
distributed, with mean and standard deviation 1. 
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Fig. 2. — Top: Plot of Ah=h — h re f as a function of energy e = his, with h re f being the 
unknown altitude of peak reference energy e re f emission. e re f was chosen to be 35 keV (see 
text for more details). The error bars were propagated from the fittings. Black: using 
all our events, blue: using all our events with maximum pileup error below 20%, magenta: 
using all our events with maximum pileup error above 20%. Bottom: Density structure 
(assuming fully- ionized plasma), as a function of Ah=h — h re f derived from the top plot (see 
text for details). Some data points at high densities have no error bars: this means the error 
bar was actually larger than the nominal value. We kept them on the plot because of their 
consistency. The squares represent the densities derived around 35 keV emission. 
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The data presented in Figure [2] (top) will be used in Section [3] to determine an average 
solar chromopheric/coronal density structure during flaring times. It can already be clearly 
seen that above ~50 keV, flare footpoint centroids have a vertical extension of less than 
~0.5 Mm, i.e. electrons above ~50 keV reach their stopping heights within that 0.5 Mm 
region. We finally discuss the potentially important sources of contamination to our dataset 
in the next subsection. 



2.3. Possible sources of contamination 



Presence of a thermal component: At low energies (e.g. usually <15 keV for M-flares, 
and <25 keV for X-flares), the non-thermal thick-target beam model can no longer be 
applied, due to the presence of high-temperature, high-altitude X-ray emitting loops. 



Presence of a non-thermal coronal source: (IKrucker et al.l 120081 ) These sources are 
often difficult to o bserve, due to instrum ental dynamic range, but are probably present 



most of the time (IKrucker &: Lin 



20081 1 . It is likely that such high-altitude source 
exist, slightly moving the center of gravity of our footpoint sources to higher projected 
altitudes. Assuming a coronal source altitude of 10 Mm above the footpoint (s), with 
10% of the footpoint non-thermal flux, the upward shift in altitude of the emission 



centroid cou 
spectra (e.g. 



d be ~1 Mm. As coron al source spectra are much steeper than footpoint 



Battaglia &: Benz 



20061 ). this effect is less pronounced at higher energies. 



2002|) 



Pulse pileup: At high photon fluxes, detector pulse-pileup occurs ( ISmith et al.l 
This has the effect of combining two low-energy photon into a single higher-energy 
photon. Hence, an intense thermal loop at ~18 keV will have a trace in ~36 keV 
images (in the case of attenuator state 3). A pile-up error of 10% (as given e.g. by the 
routine hsi_pileup_check.pro) under attenuator state 3, should roughly shift the center 
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of gravity of 36 keV emission by ~1 Mm to higher altitudes (assuming the thermal 
loop is about 10 Mm above the footpoints) This effect will typically occur only at 
high countrates, and is very energy- dependent, i.e. present mostly at twice the energy 
where most of the counts are, itself dependent on the spacecraft attenuator state: ~12 
keV for attenutator state 0, ~24 keV for attenutator state 1, ~36 keV for attenutator 
state 3. 



Albedo: Albedo effects rtBai fc Ramaty 



1978 



Kontar fc Jeffrey 



20101 ) can also influence 



the results: up to 20-40% of the 30-40 keV flux can be due to albedo (assuming 
isotropic beaming). Assuming most of the 30-40 keV emission comes from ~2 Mm 
above the photosphere, a ~30% component reflected from the photosphere would 
shift the centroid of 30-40 keV emission by about ~0.5 Mm downwards. This effect 
is present for all flares, is expected to be greatest around 30-40 keV, and shifts the 
centroid to lower altitudes, contrary to the pulse-pileup and coronal source effects. 
Furthermore, the amplitude of the effect of albedo is maximum near disc center, and 
minimum near the limb, and hence has a tendency to offset (in part) the effects of 
pulse pile-up and coronal sources. 



To conclude, it can be said that observed heights for energies >50 keV are trustworthy 
(within their statistical limitations), but that non-thermal emission <40 keV may be 
slightly offset (typically by up to 1-2 Mm) to higher altitudes than in reality, and with the 
magnitude of this offset varying somewhat with energy 

With these caveats in mind, we have attempted in the next section to derive an average 
density- height profile from our dataset. 
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3. Derivation of densities 

3.1. Direct derivation 

Figure |2] (bottom) displays the density structure obtained using all the flares in our 
list, and usi ng two subsets of it , with error bars. They were obtained using the method 

(120021 ). and in Appendix E 



described in 



Brown et al. 



A fully-ionized corona was assumed, but high energy electrons (producing most of the 
high-energy radiation) probably reach regions of low ionization in the lower chromosphere. 
This means that densities obtained at low altitudes (derived from the higher energies) 



bremsstrahlung 



be mu 


tiplied by ~2.8 (electron 


Brown 


1973; 


Kontar et al. 


2002) 



Figure [2] suggests that 1) emission at lowest energies are clearly contaminated by 
thermal loops, and results derived from those assuming non-thermal beam models are not 
to be trusted; 2) emission above 50 keV is emitted from regions with densities >10 13 cm -3 , 
assuming fully-ionized plasma. As we are probably in the chromospheric neutral region, the 
densities are likely ~3 times as much. 3) emission at intermediate energies (~20-40 keV) 
appear to come from regions of different densities, depending on the amount of pileup, itself 
dependent on the intensity of the non-thermal flux. From 2), one can hypothesize that >50 
keV emission always comes from the same altitudes, whether we are dealing with small 
or large flares. This in turn tells us that emission at ~35 keV for high pile-up events is 
almost 1 Mm above that of low pile- up events (top plot of Figure [2]), which is in accordance 



with the discussion in Section 12. 3[ but does not preclude the following scenario: larger 
(and longer) flares, which typically have high pileup, are expected to have had greater 
chromospheric ablation by the time we reach the peak in non-thermal flux (when our 
observation are made), making the loop denser, and the emission at intermediate energies 
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(~35 keV) higher. What we observe is an average of this effect over all our selected flares. 

Figure |5] (bottom) suggests an atmosphere with density profile with at least two 
exponential components. In the next section, we have formulated simple density models 
and attempted to fit it to the data in Figure E](top). 



3.2. Density model fitting 

We have attempted to fit our data with several density models, but will only discuss 
a two-exponential with unit step ionization (at altitude h s t ep ' fully-ionized for altitudes 
h > h step and fully neutral for h < h step ) model, and only briefly mention some of the others. 



As in iKontar et al.l (120081 ) . we have added a data point for our fittings, which actually 
helps us in determining h re f, our reference height for emission at e re /=35 keV: it is the 
well-established hydrogen density at the photosphere n(h = 0) = n =1.16xl0 17 cm -3 . 

We have used the routine mpfit.pro (found e.g. in the IDL Astronomy Library) to 
make our fittings (and tried also IDL's amoeba.pro routine, with no discernible differences). 
We have used a Monte Carlo approach to determine error bars for our fitting parameters: 
for each fitting "run" (a hundred such runs were executed), a random amount was added to 
each data point. This amount is normally distributed, with mean 0, and standard deviation 
equal to the nominal error of the data value. 

We show only the two-barometric component with a unit step ionization change model 
(Figure |3]) for the low pileup case, because we want to minimize the influence of pulse-pileup. 
Data with higher pileup and/or other more complicated models produced large y 2 results. 
We only fitted above ~20 keV. Events with low pileup errors have typically little to no 
thermal emission above this threshold. 
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Fig. 3. — Two-barometric component atmosphere fitting on low-pileup data (<15%), with 
unit step ionization variation. Top: Actual data points and fit (the horizontal "error bars" 



spher ic model (jVernazza et al 



1981 



19931 . solid purple), CICM (Ca l tech I rreference Chromosheric Model 



green), and lAschwanden et al. 



~ *-"V> 

model ( 


Fontenla et al. 


Ewell et al. 


I993J, solid 



(120021 ) (olive green) results (which assumed a power-law den- 
sity distribution and a fully-ionized atmosphere). Solid orange: FAL-P ionization level, 
dashed orange: VAL-C ionization level. 

Table [1] summarizes the result of the fittings. The h step parameter is the altitude where 
the ionization abruptly changes from 100% (corona) to 0% (chromosphere), and e s t ep is the 
minimal initial energy required for electrons to reach the neutral layer, hox is the altitude 
of the change from the low-altitude exponential component to high-altitude exponential 
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Table 1: Fitting parameters and other parameters derived from them. Double exponential 



structure, with unit step ionization, using flares with low pileup. 


Fitting parameters 


Description 


Units 


Value 


#0 


Low-altitude scale height 


km 


131±16 


hoi 


Altitude of transition between exponential components 


Mm 


1.69±0.21 


Hi 


High-altitude scale height 


Mm 


5.4±0.6 


hstep 


Altitude of unit step ionization change 


Mm 


1.3±0.2 


h 

1 b acc 


Altitude of acceleration region 


Mm 


19±32 


x 2 






0.79 


Derived 


^acc 


Density at acceleration altitude 


cm 3 


(8 ±71) x 10 8 


E-step 


Initial electron energy required to reach neutral layer 


keV 


36±6 



component. The scale heights H and Hi are trustworthy, thanks to the additional Uq 
datapoint and to the tight error bars at lower energies. The acceleration altitude h acc 
and density n acc are unsurprisingly poorly constrained. The altitude of ionization change, 
h s tep, is found to be «1.3±0.2 Mm above the photosphere, close to the 5%-10% ionization 
level in the VAL-C or FAL-P atmospheres (Figure [3]). The low-altitude scale height was 
determined to be ~131±16 km, which corresponds to a (neutral gas) temperature of 
5600±700 K, consistent with chromopheric temperatures below the transition region. The 
high-altitude scale height of the flaring atmosphere is determined to be ~5.4 Mm, which, 
assuming a fully- ionized isothermal plasma, corresponds to a temperature of ~115 kK. As 
this region of the lower corona is very dynamic (i.e. non isothermal) during flares, and this 
value is an average over more than 800 flares, we suspect this usually mid-transition region 
temperature does not have any intrinsic value. 

We would like to point out that events with higher maximum pulse pileup error 
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required another exponential component at intermediate altitudes in order to produce a 
"reasonable" fit (a l beit w ith reduced-^ 2 >4). This led to a fit very close to the CICM and 



Aschwanden et al. 



(120021 ) results at intermediate altitudes (1.5-3 Mm), and, at low and 
high altitudes, similar to the shown low-pileup case. While an appealing result, we do not 
think it is trustworthy, given the afore-mentioned issues. 



4. Conclusion and future work 

Using most of RHESSI's observations of flares with strong non-thermal HXR emission, 
we have derived an average emission energy vs. relative height profile. From this curve, it 
has been possible to derive densities, assuming a beam-like thick-target model. Furthermore, 
adding the well-known density at the surface of the photosphere as an additional data 
point, we could derive an average absolute height vs. density profile, and successfully fitted 
a two-exponential with step ionization model atmosphere to it. 



Although imaging show larger vertical spatial extent of a few Mm (IKontar et al. 



2010h . 



we have found that flare footpoint centroids above ~50 keV extent vertically, as expected 
from the thick-target beam model, and are within ~0.5 Mm of each other, in regions with 
neutral densities well above 3xl0 13 cm -3 . We have found density scale heights at low 



altitudes/high energies of 131±16 km, matching wel 



well as the 155±30 km value found by 



Kontar et al 



the VAL-C or FAL-P models (as 



2010|), and of 5.4±0.6 Mm at higher 



altitudes, for a flaring atmosphere. A host of other fitting parameters at intermediate 
heights or energies, and other derived quantities, could be obtained, but various intrumental, 
physical and/or statistical effects prevented us from obtaining reliable results for these. We 
have listed the different major sources of contamination and gave example of their possible 
contributions: thermal loops, non-thermal coronal source, pulse pileup, and albedo. The 
fitting results we presented minimized the impact of pulse-pileup, and were not affected by 
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thermal contamination. 

While properly accounting for all these contaminating effects is difficult in a statistical 
study (particularly the contribution of any coronal source), we believe they are quantifiable 
and can be compensated for in a few well-chosen events. But in keeping with the statistical 
approach used so far, we plan on investigating the use of backprojection-based imaging with 
a very limited set of fine subcollimators: it is our hope that the lower sensitivity stemming 
from the use of a smaller number of collimators will be offset by the fact that taking the 
brightest pixel (as opposed to a centroid) in spatially-resolved sources should be less prone 
to the contaminants we have discussed so far. Our goal is to better resolve the intricate 
interplay between changes in density and changes in ionization level at intermediate 
altitudes (~l-3 Mm above photosphere, corresponding to intermediate energies of ~ 25-50 
keV). 



A. Determining altitude differences from 2-D solar maps 

A point in space has coordinates (r,8,4>) in a spherical coordinate system centered on 
the Sun. Transforming into rectangular coordinates: 

x = rsin6>cos0 (Al) 
y = r sin 8 sin <fi (A2) 
z = TCOS0 (A3) 

Most solar maps use a planar coordinate system (X,Y) centered on the Sun, where X 
and Y represent angular distance from Sun center, as observed from Earth. Assuming that 
Earth is somewhere along the x axis, we have: 

X ■ al = y = r sin 9 sin </> (A4) 
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Y-d 



z = r cos 9 



(A5) 



Where d is the Sun-Earth distance. 

On a solar map, emission at energy e$ is located at position (ATj,YJ). In this Appendix, 
we use lower case letters for real quantities, and upper case letters for quantities projected 
on the Sun. Assuming that emission at different energies is along a radial from the Sun, we 
have: 



Xi ■ d = (r s + hi) sin # sin < 
Yi ■ d = (r s + hi) cos 9 



(A6) 
(A7) 



where r s is the solar radius, and hi the altitude of emission at energy £j 



Hence, between energies E\ and e 2 . 



AX-d = {X 2 -X 1 )-d=(h3-h 1 )an9mi 
X x + X 2 



X-d 



■ d 



,h 1 + h 2 . . , . 
r s H I sin 9 sin 



Ah sin 9 sin (j) 
i r s sin 9 sin > 



(A8) 
(A9) 



And 

X r s 



2 V 2 

, or Ah = f B ^=-. The ratio 4^ is obtained from no n-thermal HXR flare ma ps 



x 



and r s = 696 Mm, yielding Ah. This method is the one used by 



Matsushita et al. 



(1992) 



Alternatively (e.g. ISatoll2006h : 



AR-d = ^/{X 2 - X x f + (Y 2 - n) 2 = Ah V sin 2 9 sin 2 + cos 2 9 (AlO) 

(All) 



R-d 



X 2 + Y 2 m r s \l sin 2 9 sin 2 6 + cos 2 9 



And Ah=r s ^=i, with ^ given by HXR flare maps. 

In this paper, we used this second method to determine Ah. As it was using twice as 
much information, it was deemed the best. Practically, both yielded very similar results. 



See Figure [T] for how the slopes ^ are derived from linear fitting of the data. 
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B. Direct derivation of densities 



As noted in iBrown et all (120021 ) (see also ISaint-Hilaire et all 120091 ) , and for electrons 
having a negative injected power-law distribution and propagating towards region of higher 
densities, non-thermal emission at energy e is principally emitted by electrons near energy 
e, and spatially peaks where electrons have crossed a column density: 

-2 



N(s) 



f s e 1 
I n{s ' )ds ' W 2K 



(Bl) 



with s the distance from the acceleration region, e in keV, and if=2. 6x!0 cm keV 



(fully- ionized corona). It is important to note (as amply discussed in 



Brown et al. 



20021 ) 



that N(s) 



2K 



is an approximation, and that in truth a depen dence on 



inject ed electrons and on the density structure exists (see also 



the spectrum of the 



Xu et al. 



2008 



Prato et al. 



2009|). 



Using the previous approximation, and assuming that s is radial and using instead the 
altitude z as variable, the density n(z) is simply obtained through derivation: 

dN(z) e de 



n(z) 



dz 



Kdz 



(B2) 



(In a more general treatment, one can account for the changing ionization level with altitude 
by making K a function of z.) 

For discrete data points, where z% are the position of maximal emission at energy Si, 
and for i^j: 

Zi + Zj . (fTj 4" Si) Sj £j 1 Sj £j /T->n\ 

nizii = w — —— = (B3) 



2K z< 



Our data yields Ei and Ahi=Zi — z re f. Zij=h re f+ Ah ^ Ah:> , with h re f to be determined by 



other means. 
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C. Model fitting 

C.l. Single exponential with uniform ionization 

This simple analytical model will not be used, but is shown as a starting point for 
other models. 

Ionization level is assumed 100% throughout. The density structure is modeled using: 

n(h) = n e- h / H = n acc e^ h - h ^ H (CI) 

with no = 1.16 x 10 17 cm -3 , the well-accepted photospheric value, h is the altitude 
above photosphere, H is the scale height, and h acc and n acc the height and density of the 
acceleration region. 

The column density is: 
f h 

N(h) = - n{h)dh = n H ( e - h/H - e ~ hacc/H ) = n acc H ( e -( fc - fc — V H - l) (C2) 

" hacc 

where n acc is the density in the acceleration region. 
Using the e 2 ~ 2KN approximation, we get: 

e 2 = 2Kn acc {e^ h - h ^ ),H - l) (C3) 

and: 

k = k '" = - Jh ( 1 + st) (C4) 

And the difference in height Ah™ between the centroid of emission at energies E{ and 



Sf 



e 2 

1 + ^- 



Ahij = hi - hj = H In I : 2KH ^ I (C5) 



2KHn a 



The data yield Ahij, and Ej, which can be used to determine the fitting parameters 
H and n acc . And, finally, h re f can be determined using Eq. IC4I with E=E re f. Notice that 
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only good knowledge of H is required to determine h re f with accuracy, and is independent 
of what n acc (or, alternatively, h acc ) is. As we will see, this will prove to be very useful, as 
n acc or h acc will often turn out to have quite large errors. 



C.2. Double exponential with uniform ionization 

Ionization level is assumed 100% throughout. The density structure n(h) has h as 
variable and H , Hi, and h i as parameters. 

{ n e -h/H h<h m 
(C6) 
nQie -(h-h 01 )/ Hl ^ h>h 01 

with: 

n = 1.16 x 10 17 cm" 3 (C7) 
n 01 = n e- hoi/Ho (C8) 
n acc = n 01 e-( h -- hoi)/Hi (C9) 

(CIO) 

The column density iV(/i) requires an additional parameter, h acc , the height of the 
acceleration region. 

N(h) = - I n(h)dh (Cll) 

which can easily be integrated to: 



N(h) 



, h > h acc 

H inacc (e^ h - h ^y^ - 1) , h 01 < h < h a ipi2) 

^ H 1 n ae c(e- ( ' hoi ~ haa>)/Hl - 1) + H n 01 (e- ( - h - hai ^ Ho - 1) , h < h 



'01 



The rest of the procedure is similar to as explained in Section IC.lt except for the 
following: (a) We have two additional fitting parameters hoi an d Hi. (b) The explicit e(h) 



expression that we get is no longer easily invertible into a simple h(e) equation, as done in 
IC3I and IC4I The inversion is hence done numerically, using interpolations. 



C.3. Double exponential with unit step ionization 



We use the same equations as in the previous section, except we add another parameter, 
hst e p, the altitude at wh ich the solar atmop shere is assumed to abruptly change from fully 



ionized to fully neutral (IKontar et al. 



2002 



and citations therein). In the computation o 



N(h), the density in regions below h step is weighted by a factor 1/2.818 (IKontar et al.ll2002j ). 



C.4. Triple exponential with uniform ionization, or with unit step ionization 

change 

As in IC. 21 and IC.3t with an additional intermediate exponential component. 



D. Analytical proof that < 0, V7V and V5 > 2, with consequence that for an 
injected negative power-law (5 > 2) distribution of electrons propagating 
under the thick-target assumption to have peaks of increasing energy 
emission as it propagates, the density gradient along the beam's path must 

be positive: 



From 



Brown et al 



(12002) : 



Saint-Hilaire et al. 



( 120091 ). the emissivity per unit column 



density of such a distribution of electron is: 



1 8 1 

1 + u 2' 2 



(Dl) 



where e is the emitted photon energy, iV the column density crossed by the propagating 
electrons, 5 the negative spectral index of the injected electrons, (3 the incomplete Beta 
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function, u 



(see e.g. 



— 2lfiV' 




Brown et al. 


2002 



, a nd the constants K=2.6xl0" 18 cm 2 keV~ 2 , and A={5 - l)j±z 2 §^E 5 l 
for a complete explanation of these constants). 



(flat with iV) when 



Saint-Hilaire et al. 



2/-; 



^7 ^> 1, and is oc e 1 N 8 I 2 when ^ 1 ( see Fig ure HI t°P an d e -g- 



2009 



, , Appendix A). Straightforward derivation yields: 



u 



1/2/ 



1 + ti) 



-5/2-1/2 



1 5 1 

l + u'2'2 



(D2) 



which happens to have the shape shown in Figure H] (bottom): Notice that 4Uz <0 for 



1 cr = i 
10~ 6 r 
1 0" 7 ? 

io _a ; 

0. 

0- 24 F 
0" 25 r 

-26 

0- 27 r 
-2s ; 



1.0 10.0 

(2KN) 1/2 



1.0 10.0 
(2KN) 1/2 



100.0 



100.0 



Fig. 4. — Shape of 4^ and — (both in arbitrary units) as a function of ^2KN, the latest 
corresponding to the minimal initial electron energy required to cross column density N. 



all 5, e, and iV of interest. The mathematical proof, valid for 5 > 2, is as follows: 



d 2 i 

dN 2 



< if ^(^f; I) > u V2(i + w )-*/2-i/2 j Le . if > i ; with 



(ife f; I) M_1/2 ( 1 + «) 5/2+1/2 - Note that /3 |; |) = J 1+ " x & ' 2 - l {l - x)" 1 / 2 ^, which 



-5/2 



is > Jq +u x & I 2 1 dx, as x G [0, 1[. The latter can be integrated to |(1 + u) assuming 
5 > 2 (a requirement with most thick-target formulae, and supported observationally). 



Hence /3 (^; f ; §) > f (1 + <T 5/2 , leading to /(«, 5) > sj 1 + J. Hence /(«, 5) > 1 for all 
of interest (u G ]0, oof). QED. 



u 
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For a localized peak in emission e to exist, we must have: 

d 2 I d ( dl\ ( dn dl d 2 I\ 

d* " n dN{ n dN) =n {dNdN +n dN 2 ) (D3) 

= (D4) 

As > 0, n > 0, and < 0, this condition can only be achieved if > 0, i.e. when the 
density increases as the column density traversed increases, i,e. when our beam of electrons 
propagates towards regions of higher densities. 
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